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Abstract 

We develop techniques to determine the mass scale of invisible particles pair-produced at hadron 
colliders. We employ the constrained mass variable m,2c, which provides an event-by-event lower-bound 
to the mass scale given a mass difference. We complement this variable with a new variable rri2C,UB which 
provides an additional upper bound to the mass scale, and demonstrate its utility with a realistic case 
study of a supersymmetry model. These variables together effectively quantify the 'kink' in the function 
max rriT2 which has been proposed as a mass-determination technique for collider-produced dark matter. 
An important advantage of the mic method is that it does not rely simply on the position at the endpoint, 
but it uses the additional information contained in events which lie far from the endpoint. We found the 
mass by comparing the HERWIG generated m,2c distribution to ideal distributions for different masses. 
We find that for the case studied, with 100 fb _1 of integrated luminosity (about 400 signal events), the 
invisible particle's mass can be measured to a precision of 4.1 GeV. We conclude that this technique's 
precision and accuracy is as good as, if not better than, the best known techniques for invisible-particle 
mass-determination at hadron colliders. 

1 Introduction 

If dark matter is produced at a hadron collider, its likely signature will be missing transverse momentum. 
To help determine the underlying origin of the observed dark-matter, it is important to measure the masses 
of the new particle states. Mass determination is a key part to identifying the underlying theory which 
lies beyond the Standard Model. Newly discovered particles could be Kaluza-Klein (KK) states from extra- 
dimensions, supersymmetric partners of known states, technicolor hadrons, or something else that we have 
not anticipated. Models predict relationships between parameters: Supersymmetry relates the couplings of 
current fermions to the couplings of new bosons and the supersymmetric particle masses reflect the origin 
of supersymmetry breaking; masses of KK states tell us about the size of the extra dimensions. In general, 
mass determination of new particle states is central to discerning what lies beyond the Standard Model. 

If dark-matter particle states are discovered at the LHC, how will the masses of these new particles be 
measured? Determining the mass of the dark-matter particle (or other states which decay to dark matter) 
will be difficult because we don't expect the dark-matter particles to leave tracks or have a displaced vertex, 
and because we do not know the rest frame of the initial parton collision. 

If one has confidence about the specific model and process responsible for the observed events, then the 
magnitude of the process's cross-section constrains the mass. If one does not yet know the model, then 
model-independent techniques will be needed to determine properties of the newly discovered particles such 
as couplings, spin, phases, mixing angles, and masses. Proposed methods for finding the gluino phase |j| 
or spin rely on our ability to make model-independent kinematic mass measurements. Thus, model 
determination needs a suite of model- independent tools to provide initial constraints on the mass. 
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There has been much recent work in developing model-independent mass determination tools. Edges in 
invariant mass combinations provide information about mass differences or mass squared differences |^] . 
If the hadron collider accesses many new states, then we may be able to determine the masses by combining 
the relationships provided by many different invariant-mass edges of a cascade decay || ||, 0, There 
is also a series of approaches called Mass Shell Techniques (MST)Q where one uses an assumption about 
the topology and on-shell conditions to solve for the unknown masses. One MST variant assumes a long 
symmetric cascade decay chain and counts which multiplets of masses have solution to the most events [|l0| , 
another assumes the masses in two events must be equal [ fLl) , and another hybrid combines a MST with the 
information from the many edges in cascade decays [|l2) 0. There has been extensive work using the so-called 
'stransverse mass' variable, tot2, to determine the mass difference between parent particle's and a dark- 
matter candidate particle's mass given an assumed mass for the dark-matter candidate (Refs. (if, fl4||have 
more than 45 citations). If the final decay to the lightest state involves a three-body decay^j as shown in 
Fig. [j], then a 'kink' in the max mx2 fl^l B ; will occur at a position which indicates the invisible particle 
mass, as described by Cho, Choi Kim, Park (CCKP) 0, [if] and corroborated in Refs. flf, p|. The 'kink' 
in max rriT2 can be quantified by the constrained mass variable mic that we introduced in a previous letter 
fjpf . Even in models where new invisible particles are nearly massless (like the gravitino studied in 
one would rather not just assume the mass of the lightest supersymmetric particle (LSP), which is needed 
as an input to the traditional tot2 analysis, without measuring it in some model independent way. Which 
approach turns out to be best is likely to depend on what scenario nature hands us, since the various 
techniques involve different assumptions. Having different approaches also offers the advantage of providing 
a system of redundant checks. 

In Ref . |2(| , we introduced the 777,2c kinematic variable which gives an event-by-event lower bound on the 
dark-matter particle's absolute mass given the mass difference between the dark matter candidate and its 
parent. In this paper, we introduce a complementary variable rri2c,UB which gives an event-by-event upper 
bound on the same absolute mass. Our study shows that this technique rivals other invisible-particle mass 
determination techniques in precision and accuracy. 

In this paper, we provides a demonstration of the use of the variable m,2c and rri2c,UB in LHC conditions. 
The variables rri2c{M-) and rri2C,UB(M-) give an event-by-event lower-bound and upper-bound respectively 
on the mass of Y assuming the topology in Fig. |l| and the mass difference M_ = my — tun. To get the mass 
difference, we use events where Y decays into N and two visible states via a three-body decay in which we 
can easily determine the mass difference from the end point of the visible-states invariant-mass distribution, 
m\ 2 - One might also conceive of a situation with 7712c supplementing an alternative technique that gives a 
tight constraint on the mass difference but may have multiple solutions or a weaker constraint on the mass 
scale ||[^2]]. Given this mass difference and enough statistics, 777,2c 's endpoint gives the mass of Y. However 
the main advantage of the 777.2c method is that it does not rely simply on the position at the endpoint but 
it uses the additional information contained in events which lie far from the endpoint. As a result it gives 
a mass determination using significantly fewer events and is less sensitive to energy resolution and other 
errors. 

To illustrate the method, in this paper we study in detail the performance of the 7772c constrained mass 
variable in a specific supersymmetric model. We study events where each of the two branches have decay 
chains that end with a X2 decaying to a \° and a pair of opposite-sign same-flavor(OSSF) leptons. Thus 
the final states of interest contain four isolated leptons (made up of two OSSF pairs) and missing transverse 
momentum. Fig [l] defines the four momentum of the particle states with Y = %a> AT = x°, and the OSSF 
pairs forming the visible particles 1 — 4. Any decay products early in the decay chains of either branch are 
grouped into k which we generically refer to as upstream transverse momentum (UTM) . Nonzero k could be 
the result of initial state radiation (ISR) or decays of heavier particles further up the decay chain. Events 
with four leptons and missing transverse momentum have a very small Standard-Model background. To give 
a detailed illustration of the 777,2c methods, we have chosen to analyze the benchmark point PI from ]23| 
which corresponds to mSUGRA with 777 = 350 GeV, mi/2 = 180 GeV, tan/? = 20, sign(/i) = +, A Q = 0. 
Our SUSY particle spectrum was calculated with ISAJET |24| version 7.63. We stress that the analysis 



1 The title MST is suggested in Ref. 

2 The m,2c variable is a simple example of such a hybrid technique. 

3 The presence of a > 3-body decay is a sufficient but not necessary condition. Two-body decays can also display kinks 

|15L tlq] provided the decaying particles have sufficiently large transverse boosts. 
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Figure 1: We assume the two decay chains share a common end-state given in this diagram. All previous 
decay products are grouped into the upstream transverse momentum, k. 

technique employed applies generically to models involving decays to a massive particle state that leaves the 
detector unnoticed. 

A powerful feature of the mic distribution is that, with some mild assumptions, the shape away from 
the cndpoint can be entirely determined from the unknown mass scale and quantities that are measured. 
The ideal shape fit against early data therefore provides an early mass estimate for the invisible particle. 
This study is meant to be a guide on how to overcome difficulties in establishing and fitting the shape: 
difficulties from combinatoric issues, from differing energy resolutions for the leptons, hadrons, and missing 
transverse momentum, from backgrounds, and from large upstream transverse momentum (UTM) As we 
shall discuss, UTM actually provides surprising benefits. 

The paper is structured as follows: In Section ^| we review mic and introduce the new observation that, 
in addition to an event-by-event lower bound on my, large recoil against UTM enables one also to obtain 
an event-by-event upper bound on my. We call this quantity m2c,UB- Section [| describes the modeling and 
simulation employed. Section |] discusses the implications of several effects on the shape of the distribution 
including the m\i (in our case mu) distribution, the UTM distribution, the backgrounds, combinatorics, 
energy resolution, and missing transverse momentum cuts. In Section [5J we put these factors together and 
estimate the performance. We conclude in Section || with a discussion about the performance in comparison 
to previous work. 

2 Upper Bounds on my from Recoil against Upstream Transverse 
Momentum 

We will now review the definition of mic as providing an event-by-event lower bound on my . In generalizing 
this framework, we find a new result that one can also obtain an upper bound on the mass my when the 
two parent particles Y recoil against some large upstream transverse momentum kx- 

2.1 Review of the Lower Bound on my 

Fig |l| gives the relevant topology and the momentum assignments. The visible particles 1 and 2 and invisible 
particle N are labeled with with momentum ct\ and (which we group into a = a.\ +0.2) and p, respectively 
j3 = fix + /3a and q in the other branch. We assume that the parent particle Y is the same in both branches 
so (p + a) 2 = (q + f3) 2 . Any earlier decay products of either branch are grouped into the upstream transverse 
momentum (UTM) 4- vector momentum, k. 

4 Our references to UTM correspond to the Significant Transverse Momentum (SPT), pair production category in Jlrj ] where 
SPT indicates that the relevant pair of parent particles can be seen as recoiling against a significant transverse momentum. 
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In previous work with events of this topology, we pcj showed how to find an event-by-event lower bound 
on the true mass of and my. We assume that the mass difference M_ = my — mjv can be accurately 
measured from the invariant mass edges max mia or max 77134. For each event, the variable 7712c is the 
minimum value of the mass of Y (the second lightest state) after minimizing over the unknown division of 
the missing transverse momentum f T between the two dark-matter particles iV: 

af (1) 

(2) 
(3) 
(4) 
(5) 

There are eight unknowns corresponding to the four momentum components of p and q and five equations of 
constraint. Because the true p and q must be within the domain over which we are minimizing, the function 
is guaranteed to be less than or equal to my. 

One way of calculating 7712c f° r an event is to use my^Cxiv) Jl3|, which provides a lower bound on 
the mass of Y for an assumed mass xn oiN. The true mass of Y lies along the line xy {xn) = M-+XN where 
we use xy to denote the possible masses of Y and to distinguish it from the true mass of Y denoted with my. 
Thus we can see that for xn to be compatible with an event, we must have 7777-2 (xw) < Xy(xn) = Xn + M_. 

For a given event, if one assumes a mass xn for N, and if the inequality mT2(Xiv) < Xn + M_ is 
satisfied, then there is no contradiction, and the event is compatible with this value of xn- If however, 
itit2(xn) > M- + Xn, then we have a contradiction, and the event excludes this value xn as a viable mass 
of N. Using this observation, m 2 c can be found for each event by seeking an intersection between itit2(xn) 
and xn + M- p0| ■ The lower bound on my is given by my > Af_ + x°n where x°n IS a zer0 °f 

9(Xn) = m T2 (xN) -Xn - M- 

with g'(x o N )<0- (6) 

In the case k = 0, the extreme events analyzed in CCKP |0 demonstrate that g(x) will only have one 
positive zero or no positive zeros, and the slope at a zero will always be negative. For no positive zeros, 
the lower bound is the trivial lower bound given by M_. Note that a lower bound on the value of my 
corresponds to a lower bound on the value of m^r. The Appendix in Ref. |2C|] shows that at the zeros of 
9(Xn) which satisfy Eq(pj), the momenta satisfy Eqns(E 



2.2 A New Upper Bound on my 

If there is large upstream transverse momentum (UTM) {kr ~ M_) against which the system recoils, then 
we find a new result. Using the method to calculate 7712c gives one the immediate ability to see that my 
can also have an upper bound when requiring Eqns(||-|^). This follows because for large UTM the function 
9(Xn) may have two zeros[] which provides both an upper and a lower bound for my from a single event. 
We have also found regions of parameter space where g(x) has a single zero but 5'(xat) > corresponding 
to an upper bound on the true mass of m^r ( and my) and only the trivial lower bound of > 0. 

We can can obtain some insight into the cases in which events with large UTM provides upper bounds 
on the mass by studying a class of extreme event with two hard jets, j a and jp against which Y recoils 
(k = j a + jp). We will describe this extreme event and solve for the regions of parameter space for which 
one can analytically see the intersection points giving a lower bound and/or an upper bound. The event is 
extremal in that 771t2(xa0j which gives a lower bound on my, actually gives the true value of my when one 
selects xn equal to the true mass mjy. 

The ideal event we consider is where a heavier state G is pair-produced on shell at threshold. For 
simplicity we assume the lab-frame is the collision frame. Assume that the Gs, initially at rest, decay into 

5 There may be regions in parameter space where function g(x) has more than two zeros, but we have not encountered such 
cases in our simulations. 
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visible massless jets j a , jp and the two Y states with the decay product momenta a+p and (3 + q. Both jets 
have their momenta in the same transverse plane along the negative i-axis, and both V's momentum are 
directed along the z-axis. Finally, in the rest frame of the two Ys, both decay such that the decay products 
visible states have their momentum a and /? along the i-axis and both invisible massive states N have their 
two momenta along the negative i-axis. In the lab frame, the four-vectors are given by 

" {1,-1,0,0} (7) 

a-f - {1,1, 0,0} (8) 

For the event given by Eqns(^-^), we can exactly calculate itit2{xn)'- 

2x 2 N m Y + (mjf-mfy) (m 2 N m G -rn Y ) + (m Y -rn 2 N ) J A m%x 2 N m y + (m Y - m 2 N m 2 G ) 2 
m T2 (xN) = 2< ■ 

(10) 

This is found by calculating the transverse mass for each branch while assuming xn to be the mass of N. 
The value of p x is chosen so that the transverse masses of the two branches are equal. Substituting this 
value back into the transverse mass of either branch gives mx2(xjv)- 

Fig H shows g(xN), given in Eq(||), for several choices of m G for the process described by Eqs(fjj]|) 
with M- = 53 GeV and m^r = 67.4 GeV. Because G is the parent of Y, we must have mg > my. If 
m Y < m G < 2m Y /(mN + m Y ), then itit2(xn < it^n) is larger than xn + M- up until their point of 
intersection at xn = tun- I n this case their point of intersection provides a lower bound as illustrated by 
the dotted line in Fig. [2] for the case with rriQ = 150 GeV. For 2m Y /(mN + my) < m G < y/m Y /mN there 
are two solutions 

XN,Min = m N (11) 

(m^ — my) (—2m Y + mNm G mY + m 2 N m G ^j 

XiV,Max = 7 —, -Z v T7 j -. —Z TT (12) 

(mjvmG + \ m G — 2m Y )m Y )(m]^m G + m Y (m G + 2m Y )) 



When me = y/m Y /iriN, the function g(xN) has only one zero with the lower bound equalling the upper 

bound at m n . The solid line in Fig. || shows this case. Between \J m Y /vtln < m G < \/6m Y / \J (mjv + fny ){2mj^ + m Y ) 

we again have two solutions but this time with 

(m N - m Y ) {-2m Y + m N m 2 G m Y + m 2 N m 2 G ) 

XAT,Min — 7— — — —7 ^ > T7 j -, —Z 77 

{uinitig + (iTiG — 2mY)my)(mNmG + my\m,G + 2my)) 

XAT,Max = m N . (14) 

The dashed line in Fig. || shows this case with tug — 170 GeV. For greater than this, we have x_/v,Max = 
m N and XAf.Min = 0. 

This example illustrates how m-ic can provide both a lower-bound and an upper-bound on the true mass 
for those events with large UTM. The upper-bound distribution provides extra information that can also be 
used to improve early mass determination, and in what follows we will refer to the upper bound as m,2C,UB- 
We now move on to discuss modeling and simulation of this new observation. 



3 Modeling and Simulation 

As a specific example of the application of the m 2 c method, we have chosen a supersymmetry model 
mSUGRA, m = 350 GeV, m 1/2 = 180 GeV, tan/3 = 20, sign(^i) = +, A a = |. The spectrum used in 

6 This was model PI from fea| which wc also used in [M. 
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Figure 2: Shows g(xN) for the extreme event in Eq(0-^|) with M_ = 53 GeV and mjv = 67.4 GeV. The 
dotted line has = 150 GeV and shows an event providing a lower bound on my- The dashed line 
rriQ = 170 GeV and shows an event with both a lower bound and an upper-bound on my. The the solid 
line shows vtiq = \J my JmN where the lower bound equals the upper bound. 



the simulation has to^° = 67.4 GeV and m^o = 120.0 GeV. We have employed two simulation packages. 
One is a Mathematica code that creates the 'ideal' distributions based only on very simple assumptions and 
input data. The second is HERWIG |25|, ^6[ |27j which simulates events based on a SUSY spectrum, MSSM 
cross sections, decay chains, and appropriate parton distribution functions. If the simple Mathematica 
simulator predicts 'ideal' shapes that agree with HERWIG generator, then one has reason to believe that the 
all the relevant factors relating to the shape are identified in the simple Mathematica simulation. This is an 
important check in validating the benefits of fitting the m^c and rri2C,UB distribution shape as a method to 
measure the mass of new invisible particles produced at hadron colliders. 

3.1 Generation of "Ideal" Distributions 

Our 'ideal' distributions are produced from a home-grown Monte Carlo event generator written in Mathemat- 
ica. This generater serves to ensure that we understand the origin of the distribution shape. It also ensures 
that we have control over measuring the parameters needed to determine the mass without knowing the full 
model, coupling coefficients, or parton distribution functions. We also use this simulation to determine on 
what properties the ideal distributions depends. 

The simulator is used to create events satisfying the topology shown in Fig [l] for a set of specified masses. 
We take as given the previously measured mass difference — m xi = 52.6 GeV, which we use in all our 
simulations. We neglect finite widths of the particle states as most are in the sub GeV range for the model 
we are considering. We neglect spin correlations between the two branches. We perform the simulations 
in the center-of-mass frame because m^c and rn,2c,UB are transverse observables and are invariant under 
longitudinal boosts. The collision energy y/s is distributed according to normalized distribution 



Is — 4 Tfl?~ 

p(^ = 12m| £ V g2 * 2 (15) 

unless otherwise specified. The X2 i s produced with a uniform angular distribution, and all subsequent 
decays have uniform angular distribution in the rest frame of the parent. The UTM is simulated by making 
kr equal to the UTM with k 2 = (100 GeV) 2 (unless otherwise specified), and boosting the other four-vectors 
of the event such that the total transverse momentum is zero. As we will show, these simple assumptions 
capture the important elements of the process. Being relatively model independent, they provide a means 
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Figure 3: The mic and rri2C,UB distributions of HERWIG events before smearing (to simulate detector reso- 
lution) is applied. The distributions' end-points show sa 120 GeV. The top thick curve shows the net 
distribution, the next curve down shows the contribution of only the signal events, the bottom dashed curve 
shows the contribution of only background events. 



of determining the mass for various production mechanisms. If one were to assume detailed knowledge of 
the production process, it would possible to obtain a better mass determination by using a more complete 
simulation like HERWIG to provide the 'ideal' distributions against which one compares with the data. Here we 
concentrate on the more model independent simulation to demonstrate that it predicts the m^c and m,ic,UB 
distributions well-enough to perform the mass determination that we demonstrate in this case-study. 



3.2 HERWIG "Data" 

In order to obtain a more realistic estimate of the problems associated with collision data, we generate samples 
of unweighted inclusive supersymmetric particle pair production, using the HERWIG Monte Carlo program 
with LHC beam conditions. These samples produce a more realistic simulation of the event structure that 
would be obtained for the supersymmetric model studied here, including the (leading order) cross sections 
and parton distributions. It includes all supersymmetric processes and so contain the relevant background 
processes as well as the particular decay chain that we wish to study. 

Charged leptons (e^ and fi ± ) produced in the decay of heavy objects (SUSY particles and W and Z 
bosons) were selected for futher study provided they satisfied basic selection criteria on transverse momentum 
(pr > 10 GeV) and pseudorapdity (|7y| < 2.5). Leptons coming from hadron decays are usually contained 
within hadronic jets and so can be experimentally rejected with high efficiency using energy or track isolation 
criteria. This latter category of leptons was therefore not used in this study. The acceptance criterion used 
for the had ronic final state was \r)\ < 5. The detector energy resolution functions used are described in 



Section 4.2 



4 Factors for Successful Shape Fitting 

There are several factors that control or affect the shape of the mic and m2C,UB distributions. We divide 
the factors into those that affect the in-principle distribution and the factors that affect the observation of 
the distribution by the detector like energy resolution and selection cuts. 

The in-principle distribution of these events is influenced by the presence or absence of spin-correlations 
between the branches, the the ran distribution of the visible particles, any significant upstream transverse 
momentum (UTM) against which the system is recoiling {e.g. gluinos or squarks decaying further up the 
decay chain), and background coming from other new-physics processes or the Standard Model. As all these 
processes effectively occur at the interaction vertex, there are some combinatoric ambiguities. These are the 
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Figure 4: We show the mic and m%c,XJB ideal distributions for five choices of assuming the HERWIG 
generated mu and UTM distributions. 



factors that influence the in-principle distribution of events that impinges on the particle detector. 

The actual distribution recorded by the detector will depend on further factors. Some factors we are able 
to regulate - for example cuts on the missing transverse momentum. Other factors depend on how well we 
understand the detector's operation - such as the energy resolution and particle identification. 

Where the effect of such factors is significant, for example for the mi2, kr, and background distributions, 
our approach has been to model their effect on the ideal distributions by using appropriate information from 
the 'data', much as one would do in a real LHC experiment. For the present our 'data' are provided by 
HERWIG, rather than LHC events, but the principle is the same. 

4.1 Factors Affecting the In-principle Distribution 
• Mass Difference and Mass Scale 

The end-point of mic and mic,UB distributions give the mass of x%- Therefore the mass scale, m%a, is a 
dominant factor in the shape of the 'ideal' distribution. This is the reason we can use these distributions to 
determine the mass scale. Fig |] shows the iri2c and rri2C,UB distributions for five choices of assuming 
the HERWIG generated mu and UTM distributions. 

How does the shape change with mass scale? The shape is typically sharply peaked at m^c = M- 
followed by a tail that ends at the mass of . The peak at M_ is due to events that are compatible with 
m^o = 0. We say these events give the trivial constraint. Because we bin the data, the height of the first 
bin depends on the bin size. As M+/M- = (rn%° + w^)/(m^» — m^o) becomes larger, then the non-trivial 
events are distributed over a wider range and the endpoint becomes less clear. In general if all other things 
are equal, the larger the mass, the more events in the first bin and a longer flatter tail. 

The distribution also depends on the mass difference M_ which we assume has been determined. We 
expect that experimentally one should be able to read off the mass difference from the m« kinematic end- 
point with very high precision. Gjelsten, Miller, and Osland estimate this edge can be measured to better 
than 0.08 GeV || using many different channels that lead to the same edge, and after modeling energy 
resolution and background. 

Errors in the mass determination propagated from the error in the mass difference in the limit of kx = 
are given approximately by 

Smyo = 1 + Sm^o = l + + 16) 

* 2 2 V M ) 2 V M ) 
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Figure 5: Dependence of m-ic distribution on the mu distribution with k = 0. Left: The mu distributions. 
Right: The corresponding m,2c distributions. The solid curves show the case where the mu distribution 
when the three-body decay is dominated by the Z boson channel, and the dashed curves show the case where 
the mu distribution is taken directly from the HERWIG simulation. 

where SM^ is the error in the determination of the mass difference M_. An error in A/_ will lead to an 7712c 
distribution with a shape and endpoint above or below the true mass in the direction indicated by Eq(|l6|) . 

To isolate this source of error from the uncertainty in the fit, we assume that the mass difference is known 
exactly in our stated results. In our case an uncertainty of <5M_ = 0.08 GeV would lead to an additional 
Sm^o = ±0.5 GeV to be added in quadrature to the error from fitting. 

• Spin Correlations 

There are no spin correlation effects relevant to m2C if the \2 pairs are directly produced and the \2 three- 
body decay is dominated by the Z boson (e.g the sleptons are much heavier that x<D pQ l- There are also no 
spin correlations if the X2 parents are part of a longer decay chain which involves a scalar at some stage. In 
the model we are studying, there are enough vertices between the two X2 decays that any correlation is very 
likely washed out, and we can treat their decays as uncorrelated. In the simple Mathematica simulations, 
we have assumed no spin dependence in the production of the hypothetical ideal distribution. 

• Input m\2 Distributions 

The mu distribution affects the m2C distribution. Fig |^ shows two mu distributions and the corresponding 
m-2c distributions with kr = (no UTM). The solid lines show the case where the three-body decay from 
X2 to X° is completely dominated by a Z boson. The dashed line shows the case where the mu distribution 
is extracted from the 'realistic' HERWIG simulation. We can see that the m2c distribution is affected most 
strongly in the first several non-zero bins. If we were to determine the mass only from the shape of these first 
several bins using only the Z contribution for the mu difference, we would estimate of the mass to be about 
4 GeV below the true mass. This can be understood because the the shape change of the mu distribution 
effectively took events out of the first bin and spread them over the larger bins simulating the effect of a 
smaller mass. 

• Input Upstream Transverse Momentum Distribution 

As we discussed in Section ||, if there is a large upstream transverse momentum (UTM) against which the 
two Xj's recoil, then we have both an upper and lower bound on the mass scale. The left frame of Fig. ^ shows 
the UTM distribution observed in the 'realistic' HERWIG data. The right frame of Fig. |6| shows the m2c an d 
m 2 c,UB distributions for fixed UTM [k T ) of 0, 75, 175, 275, 375, and 575 GeV all with k 2 = (100 GeV) 2 . 
As we discuss under the next bullet, we also find the distribution is not sensitive to the value of k 2 . For 
kr > 275 GeV, these curves begin to approach a common shape. These are ideal m 2 c upper and lower 
bound distributions where m^y = 70 GeV and my- = 123 GeV. Notice that there is no upper-bound curve 
for the case with zero kr UTM. The UTM makes the distribution have a sharper endpoint and thereby make 
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Ideal Distributions forM y =123 GeV, M_=53 GeV 
UTM from GeV to 575 GeV 




UTM (GeV) M 2 c and M 2 c ,w (GeV) 

Figure 6: Left: The UTM distribution observed in the HERWIG simulation. Right: Ideal m^c upper 
bound and lower bound distribution for a range of upstream transverse momentum (UTM) values (fcr = 
0, 75, 175, 275, 375, 575 GeV) where m N = 70 GeV and m Y = 123 GeV. 



the mass easier to determine. This is equivalent to having a sharper kink in max tot2 in the presence of large 
UTM |§. 

How do we determine kx from the data? Because we demand exactly four leptons (two OSSF pairs), 
we assume all other activity, basically the hadronic activity, in the detector is UTM. The shape used in 
the 'ideal' distribution is a superposition of the different fixed UTM distributions, shown on right frame of 
Fig. ||, weighted by the observed UTM distribution, shown on the left frame of Fig. ^. Equivalently, we 
obtain the ideal distribution by selecting kx in the Mathematica Monte Carlo according to the observed 
UTM distribution. 

• Shape Largely Independent of Parton Distributions and Collision Energy 

In the limit where there is no UTM, then vn^c is invariant under back-to-back boosts of the parent 
particles; therefore, m,2C is also invariant to changes in the parton distribution functions. 

How much of this invariance survives in the presence of large UTM? The answer is that it remains 
largely independent of the parton collision energy and largely independent of the mass k 2 as shown in Fig [?] 
numerically. On the left frame, we show three distributions and in the right frame their difference with 2 a 
error bars calculated from 15000 events. The first distribution assumes kx — 175 GeV, k 2 = (100 GeV) 2 , 
\fs distributed via [l5|. The second distribution assumes kx — 175 GeV, k 2 = (2000 GeV) 2 , \fs distributed 
via Eq(|l5|). The third distribution assumes kx = 175 GeV, k 2 = (100 GeV) 2 , and a fixed collision energy of 
y/s = 549 GeV. 

• Backgrounds 

Backgrounds affect the shape and, if not corrected for, could provide a systematic error in the estimated 
mass. In Section || we will see that the position of the minimum x 2 in a fit to is barely affected by the 
background. The main effect of the background is to shift the parabola up, giving a worse fit. To improve 
the fit, we may be able to estimate the r?i2C and ^icjjb distribution and magnitude of the background from 
the data itself. We first discuss the sources of background, and then we describe a generic technique using a 
Dalitz-like wedgebox analysis to estimate a background model which gives approximately the correct shape 
and magnitude of the background. 

One reason we study the four-lepton with missing transverse momentum channel is because of the very-low 
Standard-Model background A previous study Q estimates about 120 Standard-Model four-lepton 

events (two OSSF pairs) for 100 fb" 1 with a f T > 20 GeV cut. They suggest that we can further reduce the 
Standard-Model background by requiring several hadronic jets. Because we expect very-little direct x% P a i r 
production, this would have very little effect on the number of signal events. Also, because Zs are a part 
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Distribution Differences 




40 60 80 100 120 140 160 40 60 80 100 120 140 160 

M 1 c and M 2 c ,ub (GeV) M 2 c and M 2 c,ub (Ge V) 

Figure 7: Shows that even with large UTM, the distribution is independent of k 2 and the parton collision 
energy. Shown are three distributions and their difference calculated from 15000 events. (1) kr — 175 GeV, 
k 2 = (100 GeV) 2 , v/i distributed via Eq@. (2) k T = 175 GeV, k 2 = (2000 GeV) 2 , ^ distributed via 
Eq@. (3) k T = 175 GeV, k 2 = (100 GeV) 2 , ^ = 549 GeV. 



Signal (•) and Background(n) 




20 40 60 80 100 120 140 

m ee (GeV) 



Figure 8: The invariant mass of the OSSF leptons from both branches of forming a Dalitz-like wedgebox 
analysis. The events outside the ma < 53 GeV signal rectangle provide control samples from which we 
estimate the background shape and magnitude. The dark events are signal, the lighter events are background. 
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of the intermediate states of these background processes, very few of these events will have mu significantly 
different from mz- 

What is the source of these Standard-Model backgrounds? About 60% is from Z-pair production events 
with no invisible decay products, in which the missing transverse momentum can only arise from experimental 
particle identification and resolution errors. This implies that a slightly stronger f T cut could further 
eliminate this background. Another 40% are due to t, i, Z production. Not explicitly discussed in their 
study but representing another possible souce of backgrounds are events containing heavy baryons which 
decay leptonically. If we assume 6-quark hadrons decay to isolated leptons with a branching ratio of 0.01, 
then LHC t,t production will lead to about 10 events passing these cuts for 100 fb _1 where both OSSF 
leptons pairs have mu < mz- 

Tau decays also provide a background for our specific process of interest. The process X2 ~ > T+T ~Xi 
will be misidentified as e~e + or about 3% of the time. Because the r decays introduce new sources 

of missing transverse momentum (V T ), these events will distort the m-ic calculation. This suggests that the 
dominant background to the x^j X2 ^ + fr+ hadrons will be from other SUSY processes. 

We now create a crude background model from which we estimate the magnitude and distribution of 
the background using the 'true' HERWIG data as a guide. We follow the suggestion of Ref. |||] and 
use a wedgebox analysis plotting the invariant mass m ee against to ma1 to supplement our knowledge of the 
background events mixed in with our signal events. This wedgebox analysis, seen in Fig. |^ for our HERWIG 
simulation, shows patterns that tell about other SUSY states present. The presence of the strips along 
91 GeV indicate that particle states are being created that decay to two leptons via an on-shcll Z. The 
observation that the intensity changes above and below = 53 GeV shows that many of the states 

produced have one branch that decays via a X2 an d the other branch decays via an on-shell Z. The lack of 
events immediately above and to the right of the (53 GeV, 53 GeV) coordinate but below and to the left of 
(91 GeV, 91 GeV) coordinate suggest that symmetric process are not responsible for this background. 

We also see the density of events in the block above 53 GeV and to the right of 53 GeV suggest a cascade 
decay with an endpoint near enough to 91 GeV that it is not distinguishable from mz- Following this line of 
thinking, we model the background with a guess of an asymmetric set of events where one branch has new 
states G, X and N with masses such that the mu endpoint is 

max m?,(odd branch) = ( m G ~ m *X m * ~ m *) = ( 85 G cV) 2 (17) 

m x 

and the other branch is our X2 decay. The masses one chooses to satisfy this edge did not prove important 
so long as the mass differences were reasonably sized; we tried several different mass triplets ending with the 
LSP, and all gave similar answers. 

We now describe the background model used in our fits. One branch starts with a massive state with 
tog = 160 GeV which decays to a lepton and a new state mx = 120 GeV which in-turn decays to to a 
lepton and the LSP. The second branch has our signal decay with the X2 decaying to x° an d two leptons via 
a three-body decay. We added UTM consistent with that observed in the events. 

By matching the number of events seen outside the mu < 53 GeV region, we estimate the number of 
the events within the signal cuts that are due to backgrounds. We estimate 0.33 of the events with both 
OSSF pairs satisfying mu < 53 GeV are background events. The model also gives a reasonable distribution 
for these events. Inspecting the actual HERWIG results showed that actual fraction of background events was 
0.4. If we let the fraction be free and minimize the x 2 with respect to the background fraction, we found a 
minimum at 0.3. 

Our background model is simplistic and does not represent the actual processes, but it does a good job of 
accounting for the magnitude and the shape of the background mixed into our signal distribution. Most of the 
HERWIG background events came from W and charginos which introduce extra sources of missing transverse 
momentum. Never the less, the shape fit very accurately and the performance is discussed in Section |[ It 
is encouraging that our estimate of the background shape and magnitude is relatively insensitive to details 
of the full spectrum. Even ignoring the background, as we will see in Section ^, still leads to a minimum x 2 
at the correct mass. 

• Combinatoric Ambiguities 
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If we assume that the full cascade effectively occurs at the primary vertex (no displaced vertices), then 
the combinatoric question is a property of the ideal distribution produced in the collisions. There are no 
combinatoric issues if the two opposite-sign same-flavor lepton pairs are each different flavors. However if all 
four leptons are the same flavor, we have found that we can still identify unique branch assignments 90% of 
the time. The unique identification comes from the observation that both pairs must have an invariant mass 
mu less than the value of the max mu edge. In 90% of the events, there is only one combination that satisfies 
this requirement. This allows one to use 95% of the four lepton events without ambiguity. The first 50% are 
identified from the two OSSF pairs being of different flavors and 90% of the remaining can be identified by 
requiring both pairs satisfy mu < max mu on which branch. The events which remain ambiguous have two 
possible assignments, both of which are included with a weight of 0.5 in the distribution. 

4.2 Factors Affecting Distribution Recorded by the Detector 

As just described, the 'ideal' in-principle distribution is created from the observed mu distribution and the 
observed UTM distribution. Wc include combinatoric effects from events with four leptons of like flavors. 
Last, we can estimate the magnitude of background events and their m 2 c and m 2 c.uB shape. We now 
modify the in-principle distribution to simulate the effects of the particle detector to form our final 'ideal' 
distribution that includes all anticipated effects. The two main effects on the m 2 c and m2c,UB distributions 
are the energy resolution and the f T cuts. 

• Shape Dependence on Energy Resolution 

Energy resolution causes the m 2 c and m 2 cjjB distributions to be smeared. Here we assume the angular 
resolution is negligible. For both the Mathematica Monte Carlo model and the HERWIG events we simulate 
the detector's energy resolution by scaling the four vectors for electrons, muons, and hadrons by 

SE e 
E e 
SEy 

E„ 
SE H 

E h 

respectively [^0|||3l|. A more detailed detector simulation is of course possible, but since we do not know 
the true behavior of any LHC detector until the device begins taking data, a more sophisticated treatment 
would be of limited value here. In practice the dependence of the ideal distribution shapes on the missing 
transverse momentum resolution should reflect the actual estimated uncertainty of the missing transverse 
momentum of the observed events. 

Smearing of the distributions decreases the area difference between two normalized distributions, thereby 
decreasing the precision with which one can determine the mass from a given number of signal events. This 
expanded uncertainty can be seen in Section M, 

The m 2 c calculations depend on the mass difference, the four- momenta of the four leptons, and the miss- 
ing transverse momentum. As the lepton energy resolution is very tight, the missing transverse momentum's 
energy resolution is dominated by the hadronic energy resolution. We model the energy resolution of the 
UTM as a hadronic jet. This significantly increases the uncertainty in the missing transverse momentum 
because hadrons have about five times the energy resolution error. 

In our Mathematica model, we represent the UTM as a single four-vector k, but in reality it will be the 
sum of many four-vectors. Because we apply the energy resolution smearing to k, if k is small the simple 
Mathematica model will have a smaller missing energy resolution. However, an events with almost UTM 
could have a large missing momentum energy resolution if it has a lot of hadronic jets whose transverse 
momentum mostly cancels. Fig ^ shows that most of the time we have considerable hadronic UTM, so this 
effect is a minor correction on our results. 

• Shape Dependence on Missing Transverse Momentum Cuts 



0.1 

0.03 
0.58 



0.003 
~E~ 



0.018 



fE~H~ E H 



0.007 



0.025 
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M 2c (GeV) 

Figure 9: The missing transverse momentum vs m^c values for HERWIG data. This shows that a f T > 
20 GeV cut would not affect the distribution for rciic > 65 GeV. 

A key distinguishing feature of these events is missing transverse momentum. To eliminate the large 
number of Standard Model events with four-lepton and with no ft T , we will need to cut on this parameter. 
Fig. |^ shows the HERWIG simulation's missing transverse momentum versus the m-ic- A non-trivial ni2c 
requires substantial rf> T . Small fj, of less than about 20 GeV only affects the iri2c shape below about 
65 GeV. The shape of the rri2c < 65 GeV therefore will require a higher fidelity model from which to train 
the shapes. Instead, we just choose to not fit bins with m 2 c < 65 GeV. 

All events near the end of uiic,ub distribution require significant f Tl therefore f T cuts will not affect 
the part of this distribution which we fit. The number of events with no non-trivial upper-bounds will also 
be affected by f T cuts. We only fit the rri2c,UB distribution up to about 233 GeV. 

5 Estimated Performance 

Determining the mass based on the shape of the distribution enables one to use all the events and not just 
those near the end point. We fit both upper-bound and lower-bound shapes to the data as described in 
the appendix |a|. As one expects, fitting the lower-bound shape more tightly constrains the mass from below 
and fitting the upper-bound shape more tightly constrains the mass from above. Combining the two gives 
approximately even uncertainty. We calculate ideal distributions assuming at five values 50, 60, 70, 
80, and 90 GeV. We then fit a quadratic interpolation through the points. Our uncertainties are based on 
the value where % 2 increases by 1 from its minimum of this interpolation. This uncertainty estimate agrees 
with about 2/3 of the results falling within that range after repeated runs. Our uncertainty estimates do 
not include the error propagated from the uncertainty in the mass difference (see Eq(pj|)). 

We present results for an early LHC run, about 100 fb , and for the longest likely LHC run before an 
upgrade, about 400 fb -1 . After about 100 fb -1 , we have 700 events ( about 400 signal and 300 background). 
After 400 fb -1 , we have about 2700 events (about 1600 signal and 1100 background). Only 4 events out of 
1600 are from direct pair production. Most of our signal events follow at the end of different decay chains 
starting from gluinos or squarks. The upstream decay products produce significant UTM against which the 
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Figure 10: The result of \ 2 fits to the data with differing assumptions for 100 fb -1 (left pane) and 
400 ftV 1 (right pane). The thick line with filled squares shows the final result with all cuts, resolu- 
tion error, combinatorics, and backgrounds included and estimated in the shape fitting. This gives us 
to^o = 63.2 ± 4.1 GeV with 700 events (signal or background) representing 100 fb -1 . After 400 fb -1 this 
improves to = 66.0 ± 1.8 GeV. The error-free best case gives = 67.0 ±0.9 GeV. The correct value 
is m %o = 67 A GeV. 



two %2 parent particles recoil. 

First for the ideal case. After 400 fb -1 , using only signal events and no energy resolution, the \ 2 fits 
to the predicted shapes give = 67.0 ± 0.9 GeV (filled circles in Fig. |o|). This mass determination can 
practically be read off from the endpoints seen in Fig. I the m 2 c endpoint is near 120 GeV and subtracting 
the mass differences gives =120 GeV — M_ = 67 GeV. We now explore how well we can do with fewer 
events and after incorporating the effects listed in Section ||. 

How does background affect the fit? If we ignore the existence of background in our sample, and we fit 
all the events to the signal-only shapes, then we find a poor fit shown as the empty circle curve in Fig. [l^. 
By poor fit, we mean the x 2 is substantially larger than the 72 bins being compared (36 bins from each the 
upper-bound and lower-bound distributions). Despite this worse fit, the shape fits still give a very accurate 
mass estimate: m s? = 65.4 ± 1.8 GeV after 100 ftT 1 and = 67.4 ± 0.9 GeV after 400 fb -1 . At this 
stage, we still assume perfect energy resolution and no missing transverse momentum cut. 

Next, if we create a background model as described in Section ||, we are able to improve the \ 2 fit to 
nearly 1 per bin; the mass estimate remains about the same, but the uncertainty increases by about 20%. 
We find a small systematic shift (smaller than the uncertainty) in our mass prediction as we increase the 
fraction of the shape due to the background model vs the signal model. As we increased our fraction of 
background, we found the mass estimate was shifted down from 66.5 at 0% background to 65.6 when we 
were at 60% background. The best x 2 fit occurs with 30% background; which is very close to the 33% we 
use from the estimate, but farther from the true background fraction of about 40%. With 400 fb -1 of data, 
the systematic errors are all but eliminated with the endpoint dominating the mass estimate. These fits are 
shown as the triangles with dashed lines and give m^o = 65.1 ± 2.4 GeV which after the full run becomes 
m^o = 67.3 ± 1.1 GeV. 

Including energy resolution as described in Section ^] shows a large increase in the uncertainty. The 
dashed-line with empty square markers shows the \ 2 fit when we include both a background model and the 
effect of including energy resolution. These fits are shown as the empty squares with dashed lines and give 
m^o = 63.0 ± 3.6 GeV which after the full run becomes m-^j = 66.5 ± 1.6 GeV. 

The final shape factor that we account for are the cuts associated with the missing transverse momentum. 
After we apply cuts requiring f T > 20 GeV and fit only m-ic > 65 GeV we have our final result shown by 
the thick lines with filled squares. This includes all cuts, resolution error, combinatorics and backgrounds. 
We find m^o = 63.2 ±4.1 GeV with 700 events (signal or background) representing 100 fb -1 , and after 
400 fb -1 this improves to rn^° = 66.0 ±1.8 GeV. The true mass on which the HERWIG simulation is based 
is = 67.4, so all the estimates are within about 1 a of the true mass. 
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Figure 11: HERWIG data for 100 fb~ (thick line) and the smooth ideal expectation assuming = 70 GeV 
generated by Mathematica with all resolution, background, and combinatoric effects included (thin line). 
The x 2 °f this curve to the HERWIG gives the solid-black square on the left frame of Fig. [Ic]. 

Fig [tl] shows the ideal curve expected if m^o = 70 GeV including all effects from energy resolution, 
background, combinatoric, and f T cuts. The x 2 corresponds to the solid square on the left pane of Fig. |Tc| . 

The error in mass determination obtained with limited statistics can be estimated using Poisson statistics. 
In our studies we find that, as one would expect, increasing the number of events by a factor of four, we bring 
down our error by about a factor of two. This means that one could expect ±8 GeV after about 25 fb _1 
which represents 100 signal events and 75 background events. 



6 Discussion and Conclusions 

Despite adding some of the complicating effects one would encounter with real data, we have discovered 
other factors which demonstrate one could obtain an even better precision than we originally reported in 
Ref. ^0|. There, we used only our simple Mathematica model that assumed k = and neglected most 
sources of realistic uncertainty. We assumed all the events could be modeled as being direct production 
without spin-correlations. With these simplifications, we argued the mass could be determined to ±6 GeV 
using 250 signal events. 

In this paper, we performed a case study to show that that the relevant 7712c and rri2C,UB shapes can 
be successfully determined from the mass difference M_, the ran distribution observation, the upstream 
transverse momentum (UTM) distribution observation. We included and accounted for many realistic effects: 
we modeled the large energy-resolution error of hadronic jets. We also included the effects of backgrounds, 
f T cuts, and combinatorics. Our signal and backgrounds were generated with HERWIG. We discussed how 
a Dalitz-like plot can estimate the background fraction and shape. Observed inputs were used in a simple 
model to determine ideal distribution shapes that makes no reference to parton distribution functions, cross 
sections, or other model-dependent factors that we are not likely to know early in the process. 

Despite these extra sources of uncertainty, we found a final mass determination of ±4.1 GeV with about 
400 signal events which is still better than the appropriately scaled result from Ref. [^0| . The sources of the 
mass determination improvement are twofold: (1) the prediction and fitting of upper-bound distribution, 
and (2) the sharper end-point in the presence of large UTM. Under equivalent circumstances, the sharper 
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endpoint is enough to give a factor of 2 improvement in the uncertainty over the direct production case 
assumed in [ p0[ . Fitting the upper bound tends to improve the determination by an additional factor of y/2. 
This improvement is then used to fight the large hadronic-jet energy resolution and background uncertainty. 

Mass determination using m-ic and m2c.UB applies to many other processes. We have focused on cases 
where the mass difference is given by the end-point of an mu distribution involving a three-body decay. If 
there is not a three-body decay, then the mass difference may be found by applying other mass determination 
techniques like the mass shell techniques (MST) |l(| 11, |l2j or edges in cascade decays || [?], || or tot2 at 



different stages in symmetric decay chains [ |22| . 

How does our method's performance compare to previous mass determination methods? Firstly, this 
technique is more robust than the max m,T2 'kink' because in fitting to the shape of the distribution, it does 
not rely entirely on identification of the events near kinematic boundary. One can view m-2c and rri2C,UB as 
variables that event-by-event quantify the 'kink'. Other than the 'kink' technique, the previous techniques 
surveyed in the introduction apply to cases where there is no three-body decay from which to measure the 
mass difference directly. However, each of those techniques still constrains the mass difference with great 
accuracy. The technique of ||, @, || which uses edges from cascade decays determines the LSP mass to 
±3.4 GeV with about 500 thousand events from 300 fb _1 . The approach of jll| assumes a pair of symmetric 
decay chains and assumes two events have the same structure. They reach ±2.8 GeV using 700 signal events 
after 300 fb _1 , but have a 2.5 GeV systematic bias that needs modeling to remove. By comparison, adjusting 
to 700 signal events we achieve ±2.9 GeV without a systematic bias after propagating an error of 0.08 GeV 
in the mass difference and with all discussed effects. Uncertainty calculations differ amongst groups, some 
use repeated trial with new sets of Monte Carlo data, and others use \ 2 ■ Without a direct comparison 
under like circumstance, the optimal method is not clear; but it is clear that fitting the m 2 c and vn,-xc,UB 
distributions can determine the mass of invisible particles at least as well, if not better than the other known 
methods in both accuracy and precision. 

In summary, we have developed a mass determination technique, based on the constrained transverse 
mass, which is able to determine the mass of a dark-matter particle state produced at the LHC in events 
with large missing transverse momentum. The TO2C method, which bounds the mass from below, was 
supplemented by a new distribution m2C,UB which bounds the mass from above in events with large upstream 
transverse momentum. A particular advantage of the method is that it also obtains substantial information 
from events away from the end point allowing for a significant reduction in the error. The shape of the 
distribution away from the end-point can be determined without detailed knowledge of the underlying model, 
and as such, can provide an early estimate of the mass. Once the underlying process and model generating the 
event has been identified the structure away from the end-point can be improved using, for example, HERWIG 
to produce the process dependent shape. We performed a case-study simulation under LHC conditions 
to demonstrate that mass-determination by fitting the m^c and rri2C,UB distributions survives anticipated 
complications. With this fitting procedure it is possible to get an early measurement of the mass - with just 
400 signal events in our case study we found we would determine to^° = 63.2 ±4.1. The ultimate accuracy 
obtainable by this method is = 66.0 ± 1.8 GeV. We conclude that this technique's precision is as good 
as, if not better than, the best existing techniques. 
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A Appendix A: Least squares fit 

In order to determine we perform a \ 2 fit between ideal distributions and the HERWIG data. First for 
definitions: We define Nlb as the number of m2c events in the region to be fit, and likewise Njjb is the 
number of ni2c,UB events in the region to be fit. The m,2C of the events are grouped into bins; Cj counts 
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the events in the j bin. The variable fhsimic j-, m x°) is the normalized m,2c distribution of ideal events 
expected in bin j as calculated with an assumed m x ° , the measured M_ , the observed mu distribution, the 
observed UTM distribution, and the appropriate detector simulator. We likewise define the upper bound 
distribution to be fuB{m2C,UB ji m x°)' ^ e a ^ so define the background distribution for lower-bound and 
upper-bound distributions to be fB,LB(jn 2 c j) and / b ,u B^jnic ,u b j) and the fraction of the total events we 
estimate are from background A. 

Assuming a Poisson distribution, we assign an uncertainty, <jj, to each bin j given by 

°LB,j{™ x l) = g ( N LB ((1 - X)fLB(m 2 Cj,Tn x o) + \f B ,LB(m 2 Cj)) + Cj) , (21) 

and likewise for the upper-bound distribution. The second term has been added to ensure an appropriate 
weighting of bins with very few events that does not bias the fit towards or away from this end-point. In 
bins with few counts, normal Poisson statistics does not apply[|. 
The x 2 is given by 

X 2 (m-o) = \^ ( ° 3 ~ Nlb (1 Z A )/^g( m 2Cj' m x?) ~ N LB \ fBXB(m 2 Cj,m^) \ 2 



bin j 



E 

bm j 



Cubj - N UB (1 - A) fuB(rn 2 c,UB j ,m x o) - N jjb A /b.ub ( m 2C,UB ? - ; m x° 1 ) 



2 



&UB,j 



We calculate ideal distributions for m^o = 50,60,70,80,90 GeV. We fit quadratic interplant through the 
points. The minimum X 2 ( m x°) °f ^ ne interplant is our estimate of m^°. The amount m-^o changes for an 
increase in x 2 by one gives our 1 a uncertainty, Sm^, for J32|. 
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